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H ! Abstract 

C3 , 

| Assuming that right-handed neutrinos exist and propagate in some large extra dimen- 

q ■ sions, we attempt to give a comprehensive description of the phenomenology of neu- 

■ trino oscillations. A few alternative explanations of the atmospheric neutrino anomaly 

. emerge, different from the standard — > v T or — > v stCTi \ c interpretations. Constraints 

from nucleosynthesis and supernova 1987a are discussed. The constraints from SN1987a 
\ indicates a maximum radius of any extra dimension of about 1 A. 

<n : 
o . 

o ■ 

o ; 

^ ■ 1 Introduction 

a: 

q h ' The hypothesis that Standard Model (SM) singlets propagate in some extra dimensions with relative large radii 
D , leads to striking consequences. Most notable among them, if applied to the obvious candidate, the graviton, is the 
' possible disentanglement of the Planck scale from the scale where gravity becomes strong [BJ. It seems possible, in 
^ . fact, that more surprises from the intense theoretical activity on this and related subjects have yet to come. 

In spite of this last remark, one can start asking if it will ever be possible to make any experimental observation 
related to these phenomena. Always insisting on the graviton case, this has been and is being discussed both in 
the case of particle collisions at high energy and of tests of classical gravity Q . With the relevant parameters 
at the boundary of the allowed region, such observations, although not easy, may not be impossible. 

After the graviton, the most natural candidate to propagate in some large extra dimension is the right handed 
neutrino. Interestingly enough, the smallncss of the neutrino masses, of Dirac type, could in fact be a manifestation 
of this very hypothesis 0, @, fj], 0] . The purpose of this paper is to make a first tentative exploration of the 
related experimental consequences. A specific suggestion along these lines has already been made 

The plan of the paper is as follows. In section 2 we define the framework, including the flavour aspects. The 
possible connection with gravity is briefly outlined in section 3. In section 4 wc diagonalizc the neutrino mass 
matrix for arbitrary values of the compactification radius in presence of matter effects. In section 5 wc describe 
the oscillation amplitude of an interacting neutrino with its Kaluza-Klein tower. In section 6 we show how the 
phenomenology of neutrino oscillations gets modified. Possible signatures in on-going or future experiments are 
summarized in section 7. In section 8 we show that big-bang nucleosynthesis is not manifestly inconsistent with the 
outlined phenomenology. A strong constraint appears to come from the consideration of the neutrino luminosity 
of the supernova 1987a, as discussed in section 9. Such constraint indicates an upper limit on the largest radius of 
any extra dimension close to lA, which, if taken at face value, would make impossible any observation of oscillation 
phenomena discussed in this paper. Our conclusions are summarized in section 10. 

We view this work as a contribution to the present discussion about possibilities of physics beyond the SM, 
alternative to supersymmetric unification, which we still consider as the relatively more likely description of nature 
at small distances. 
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2 The framework defined 



Effectively, the action which defines our framework involves 5-dimensional massless fcrmions ^ i (x fi ,y), the 'right- 
handed neutrinos', one per generation, interacting on our brane with the standard left-handed neutrinos in a 
way that conserves lepton number. The relevant action is 

S = J d A x dy [^iT A id A ^i] + J d 4 x [viidvi + z^Ay ^{x^ 0)h + h.c] (1) 

where A ~ {0, . . . ,4}, xpj is one of the two independent Weyl spinors that compose fyj = (ipj,ipj), ft, is a normal 
Higgs doublet in four dimensions and A is a matrix of Yukawa couplings with dimension of (mass) . As manifest 
from (Q), A can be made diagonal without loss of generality at the only price of introducing the usual unitary 
matrix V which describes flavour changes in the charged current neutrino interactions Redefining the vi in 
terms of the neutrino flavour eigenstates 



,(/) - v . 



a = e, n, t, 



(2) 



after compactification of the 5-th dimension with suitable boundary conditions, the action ([!]) becomes S 
J d 4 xJ2 3 l= i A, with 



Ci = v>iit 



E 



Eo(a c "-or")+t.c, 



(3) 



where R is the compactification radius and nij are proportional to the eigenvalues of the Yukawa matrix A. Within 
this framework, R is the only extra parameter that will be introduced in the neutrino phenomenology other than 
the 3 masses mi and the mixing matrix V . 

As shown in appendix [A|, the Dirac nature of the mass matrix in ([}]) is made explicit by an appropriate field 
redefinition, leading to the mass Lagrangian £ ma ss = J2i=i N^MiNm + h.c. with 
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(4) 



To describe the evolution of a neutrino state of energy E u , the relevant Hamiltonian is (one for any index i) 

1 
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(5) 



diagonalized by a matrix as f/W?{M{7W T = H^ iag . The elements J7g„ give the composition of the i>i state in 



terms of the mass eigenstates corresponding to the eigenvalues A„ /2E U R 2 of TiS 

Postponing the explicit calculation of A 2 and Uq u , the amplitude A a p(t) 
t a neutrino with flavour /3, born at t = with flavour a, is, using (^|), 



A aP {t) = Yl VtVfrMt) 



(t)) for finding at any time 
(6a) 



i=l,2,3 



where 



A{ Vl Vi ,t) = ]T U oT cxp[iA« 2 t/2S,i? 2 

n=0 



(6b) 



Note the formal identity of (3a) to the standard oscillation amplitude with the replacement Aj(i) — > cxp(im 2 t/2E l/ ). 

Quite clearly, the neutrino phenomenology from (Q) can be significantly different from the standard one, but 
only if l/R is not too much bigger than all the m^. Otherwise, for l/R ^> m^, all of the Kaluza-Klcin (KK) 
states tp", ip~ n with n > 1 in (||) decouple and one remains with the standard Dirac masses miViV % R only. In turn, 
if we want to influence the atmospheric or the solar neutrino anomalies, l/R should not exceed 1 eV or so, i.e. 
i?> O.l^m. 
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3 Connection with gravity 



Before going further, let us discuss the possible connection with gravity Although not necessary, we make 
the simplifying assumption that the graviton propagates in the same extra-dimensional space as the right handed 
neutrinos. Furthermore the Planck mass Mp\ is related to the fundamental scale of the theory Mf and to the volume 
V of the (5 compactified extra dimensions as Mpi = Mf(Mj V) 1 ! 2 . To maintain a simple connection with gravity, 
we are in fact lead to consider several extra dimensions of different radii Ri, with i?i>i <C R\, so that (||) is only 
an approximation where we neglect the heavier KK excitations in the extra dimensions with small radii i?i>i ||]. 
All this is a necessity if R\ > O.lyum, as explained, and, for the heaviest neutrino mass, m max ~ v/(Mj V) 1 / 2 , as it 
follows from a simple generalization of (jl]) to S extra dimensions || . In this case 



M f M f T-r /O.lcVX^ 2 , 1 

777 ~ V- 



_ .leV— ^- and v^TJr.^I^L) lo2 6- 155 - (?) 
M P1 10 3 TeV 11 \m max J e V 5 V ; 

which shows, for any 5 > 2, the asymmetry of the different radii. 



4 Diagonalization of the neutrino Hamiltonian 

Making reference for the details to appendixH, we describe in this section the diagonalization of the Hamilto- 
nian (||), extended to include matter effects JTo[ ]7 essential at least to discuss nucleosynthesis. It is 

Wetter = ^} M ^ + ^g(^,0, 0, . . .)] (8) 

where fXi/R 2 = 2E 2 (nf — 1) depends on the refraction indices 77^ for Vi (nf) and i>i (n~) in the relevant medium. 
To be precise, the factorization of flavour as in eq.s (^) is no longer exact in presence of matter effects. Since we 
will only study flavour mixing between fx and r neutrinos in media that do not distinguish \x from r flavour, we 
can ignore this complication in the present work. 

Still calling An'' 2 the eigenvalues of 2i?„i? 2 7i^ attor , they satisfy the eigenvalue equation (see appendix A) 

A 2 — /1 — ttA£ 2 cot ttA = (9) 

where £ = mR and the index i is left understood. Hereafter R is set to unity unless explicitly reintroduced. At the 
same time the matrix elements Uo n are given by 

772 = ? (1(11 

U ° n l + n 2 e + v/\l + (A n -/VA n ) 2 /r ( ' 

Both (||) and ([h]) extend the results of Q for fi = 0. Note that, for negative fx, one A 2 eigenvalue can be negative, 
so that the corresponding A is imaginary. A qualitative description of A„ and Uo n is as follows, depending on the 
sign of \x. 



fx > case. For £ — > the eigenvalues A tend to the positive integers, except for a special eigenvalue at A 2 = fx, 
whose eigenvector tends to coincide with the ^-state up to mixings of order When fx varies with the medium 
density or temperature and crosses a positive integer, a resonant MSW conversion can take place. For £ 3> 1 
the ^-state is spread over a larger number of levels, mostly those with A « ^/JH 7r£ 2 , each with a small mixing 
w l/7r 2 £ 2 and eigenvalues close to semi-integers. Matter effects are not negligible if li>(it£ 2 ) 2 and do not 
suppress oscillations when fx is large. 



fx < case. The eigenvalues are as before, except for the special eigenvalue related to fx. For fx < — £ 2 the special 
eigenvalue A 2 , becomes negative and tends to fx as /1 -C — £ 2 - In this limit, the corresponding eigenstate is an almost 
pure v. This means that, for large negative fx, there is matter suppression. By carefully expanding in £ 2 /\/— fx, for 
the component of the instate on the special cigenstate, one finds Uf )Q ~ 1 — 7r£ 2 /2^/— \x. 
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fig- 0a fig. |b 




10" 6 1(T 4 10 -2 10° 10 2 0.01 0.1 1 10 

L/E in units of R 2 % 

Figure 1: Oscillation probabilities. Fig. 0a: P vu as junction of LjE v R 2 for £ = 1/10,1/3,1,3,10 (continuous 
lines). For £ = 1,1/3 we also show the big £ approximation (dotted lines). Fig. 0&; P vv as function of £ for 
L 3> E^R 2 (continuous line) and the small £ and big £ approximations (red dotted lines). 



5 Oscillation amplitudes for £>1 



For small £, i.e. 1/R m, only few of the KK states mix with the standard neutrino and, furthermore, with small 
mixing, of order £. In this case it is not difficult to obtain the neutrino oscillation amplitudes. 

Let us consider on the contrary the large £ limit, when the standard neutrino mixes with about 7r£ 2 KK states. 
Matter effects are neglected for the time being. For large £, A„ « n and we can safely approximate the sum in (6b) 
with an integral. Using also, from (fol), U 2 n w 2/(tt 2 £ 2 + n 2 /£ 2 ) we have 
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= e z (l — erf \fz) where z — — i ( 
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Note that the oscillation amplitude, obtained in the large £ limit, depends only on the combination £ 2 /i? 
rather than on m and R independently. For the probability P vv — \A\ 2 one has at large LjE v 



P VV {L — » oo) 

A good analytic approximation for P vu at any L is 



dn U 0n = 



1 



(ii) 



m 2 R 



(12) 
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1 — erf 
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Fig. shows P vv as function of L/E^ in units of R 2 for different values of £. The comparison of the numerical 
result with the analytic approximation at £ = 1, 1/3 shows the validity of the analytic expression even at moderate 
values of £. When £ gets sizable, the disappearance of the standard neutrino into the tower of KK states becomes 
increasingly important. Matter effects on the oscillation amplitudes can be discussed along similar lines (see 
appendix A). 



6 Unconventional fits of atmospheric and solar neutrino anomalies 

In this section we look for a comprehensive description of the neutrino anomalies, both atmospheric and solar. 
The new scale 1 / R of KK excitations allows in many different ways to have three different neutrino squarcd-mass 
splittings that could, in principle, be associated with solar, atmospheric and LSND oscillations. However, as in 
the standard phenomenology based on 3 neutrinos, we have not been able to also account for the LSND result. 
A summary of the alternative possibilities that we have found is given in table 0, as we now discuss. They are 
characterized by the different ranges of 1/R. To keep things simple, we consider a hierarchical neutrino spectrum, 
with m3 > TO2 > mi . Different spectra are only possible if all £ are small. 



4 



case £3 1/R solar 



atmospheric 



A 


> 1/3 


< 3 lO-^eV 


only VO 


— * ^ r , ^KK 


w 3 10" 3 eV 


SAM ^ e -> v KK 


with £ 2 /R w 0.01 eV. 






> 3 10" 3 eV 


standard 


B 


«l/3 


w lO^eV 


standard 


i/ M -> v r ,i/ KK or 










C 


< 1/3 


^10 _1 eV 


standard 


standard 1/^ — > z^ T 



Table 1: List of possibile oscillation patterns. 



6.1 Atmospheric neutrinos 

When 1/R decreases from about 1 eV the standard picture of neutrino oscillation is progressively perturbed. In 
cq.s (@) only A 3 is modified from the conventional form. An increasing portion of the atmospheric neutrinos starts 
oscillating into their KK towers at the expenses of the conventional — > v T transition (case C in table [l]). If 
£1,2 *C £3 *C 1 approximate formulae for the transition probabilities in an intermediate L/E v range, i.e. l/Amfj ^> 
L/E v > R 2 , are 

2tt 2 7tt 4 
P Q a = l-— l^sl^I, and for a £(3 P a[j = —\V a3 Vp 3 \ 2 & (13) 

Note that these formulae do not allow to interpret 1/R as the extra mass scale useful to account for the LSND 
result. The situation in our case, with only one extra parameter R, is far more constrained than in the case where 
one adds one sterile neutrino with arbitrary mixing parameters. 

To understand what happens when £3 increases, let us consider the case when £3 gets large. As shown in fig. ^, 
for £3 > 1 and for sufficiently large L/E v the survival probability of v 3 drops below 10%. In absence of flavour 
mixing, this would certainly disagree with the SuperKamiokande (SK) results [jfl). In cq.s with "C 1 and 

neglecting matter effects A\ and A2 are accurately approximated by the standard phase Ai = cxp(tm 2 L/2i?„), 
which reduce to 1 at the L/E u relevant to the atmospheric neutrinos. On the other hand, for £3 > 1 and above 



the SK threshold for v ^ disappearance, we can take A 3 = 0. From (6a) and above the SK threshold, using the 
unitarity of the V matrix, we have therefore 

P aa ~ (1 - |y a3 | 2 ) 2 , and for a P a0 ~ \V a3 V 03 \ 2 . (14) 

Since the data suggest P MAI w 0.5 and P ee f» 1, a fit may be possible for V e3 w and \V 2 3 \ w 0.4. These 
considerations are confirmed by an explicit fit of the dataQ shown in fig. |3|a. As mentioned in the previous section, 
for large £, the oscillation amplitudes depend on the combination £ 2 / R, which becomes the only relevant parameter 
other than V^, since we set V e3 = 0. The contour plot of the % 2 in these parameters is shown in fig. |^a. The value 
of £ 2 /i? around 0.02 eV is selected by the shape of the angular distribution of the muon neutrino events. The fit 
does not fix the value of 1/R, which can vary below about 10 -2 eV since £ > 1, i.e. m 3 > 1/R. This is case A in 
table 0. 

We can now ask what happens when £3 is decreased below 1, in the intermediate region between A and C. Now 
£3 and 1/R, as well as the mixing angles, independently influence the fit. In fig. [pD we show the profile of the x 2 
again in the plane (£ 2 /R, | 1 2 ) at fixed £3 = 1/2. At this intermediate value of £3 about half of the interacting 
neutrinos oscillate into the KK states. As a consequence a fit of the atmospheric data is possible with V^ 3 w 1, as 
illustrated in fig. ||b. Unlike the case of very large £, or of the standard fit with £ = 0, shown in fig. ||c, no — ► v T 
oscillation is present in this case. 



lr The fit of SK atmospheric data is done as in |l2j ]. In particular we have used: (a) The prediction of atmospheric v 
fluxes of 113l; (b) The energy spectra of the parent atmospheric neutrinos, corresponding to the various classes of events 
measured at SK as given by the Monte Carlo simulation of the SuperKamiokande detector, available at the www address 
www.awa.tohoku.ac.jp/-etoh/atmnu/enu/index.html. (c) The latest SK data (848 days of exposure), as extracted from (d) 
The x 2 function denned in Jij] , using 30 bins of experimental data and including all various systematic uncertainties. 'Official SK 
fits' employ unpublished data with finer energy and zenith-angle subdivisions of the neutrino-induced events. If systematic errors are 
sufficiently low, these data could discriminate the non-standard L/ E v dependence of the oscillation probability from the standard one. 
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|a: £ » 1 



|b: £=1/2 




Figure 2: /if of the SK atmospheric data in (V^ 3 , £ 2 /i?). Fig. |^a: £ 3> 1 and any i?. Fig. [|&: £ = 0.5. The contour 
lines correspond to x 2 = {15,20,25}. The best standard fit in terms of — > v T gives minx 2 ~ 14. 

Back to fig. |2|b, the two minima of the \ 2 distribution are clear. We have just discussed the case as 1. The 
other minimum originates from the one encountered before at large £. This should give an idea of what happens 
at various intermediate values of £. In this case (B in tabic 0) l/R ranges around 10 2 eV. 

This completes our discussion of atmospheric neutrinos. The fit at £ = 1/2 includes earth matter effects, 
whereas the one for large £ does not because in this case it is difficult to compute neutrino propagation across few 
layers with different density. We expect that the main features of solution A will remain unchanged. 

6.2 Solar neutrinos 

Given these alternative descriptions of the atmospheric neutrino anomaly, we should ask now how the solar neutrino 
deficit can be accounted for. So far m-i and m\ have not been fixed. We only required that £2 and £1 are both 
small. This being the case, it is in fact simple to see that a standard description of solar neutrinos is possible in all 
cases with negligible interference of the KK towers. For small values of l/R there is however the possibility of a 
transition v e — > z^kk using the MSW effect, which is compatible with the solar data It requires l/R as 3 10 -3 eV 
and a mixing with the KK states determined by £2 as 0.01, or ni2 as 10~( 4H ~ 5 ) eV, so that a fit of SK atmospheric 
data requires £3^2. When the parameter \i of section^ is specified for the electron neutrino and with the solar 
density profile, the resonant MSW conversion mentioned there (// positive, small £) takes place and suppresses the 
different components of the solar v e spectrum as possibly observed by the various solar neutrino experiments. 



7 Special features of the proposed solutions 

Some alternative descriptions of the atmospheric neutrinos appear possible. The crucial point, however, would be 
to indicate precise signatures of such solutions visible in appropriate neutrino experiments. To this purpose fig.s ^ 
arc of interest. We give there, versus L/E v , the probabilities P^ and P^ T that correspond to the fits of the SK 
results shown in fig. |^. A few features of these plots might be relevant for an experimental discrimination of the 
various possibilities. 

1. The absence of a first clear dip in the F/F„-shapc of P MA1 is a characteristic of the KK fits that we have 
discussed at intermediate and big £, at clear variance with the shape of P^ in the standard — ► v r 
interpretation of the data. 

2. The non-standard transition from unoscillated to oscillated atmospheric neutrinos requires a L/E^-range 
longer than the standard one and even the one that would be produced by neutrino decay Therefore, 
unlike what happens in the standard case, a good fit of atmospheric data significantly constrains the outcome 
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Figure 3: Best fits of the SK zenith-angle distribution of sub-GeV (first column) and multi-GeV (second column) 
VaVfi events and of upward- going muons (third column). Each plot is drawn in the (cos i9 ze nith, number of events) 
plane. The arrows on the horizontal axes denote the direction of the scattered leptons. Continuous lines denote fit 
predictions and gray bars denote the no- oscillation predictions. f„ data are plotted in blue, v e data are plotted in 
red and only in the last row, since they are the same in all rows. Crosses denote experimental data, including only 
statistical uncertainties. The systematic uncertainty in the overall number of neutrino events in each sample of 
events has been used to optimize the visual appearance of the standard fit in fig. [|c: as usual the most significant data 
are the shapes of the individual zenith-angle distributions. Fig. [|a.' fit with £ 3> 1 (v^ — ► i^kk, v t : £ 2 / R = 0.015 eV 
and V M 2 3 = 1/2). Fig. |&: fit with £ = 1/2 (v^ -> v KK : £ 2 / R = 0.004 eV and V^ 3 = I). Fig. |c: standard fit 
fa -> v T : Ato^ 3 = 3 10" 3 eV 2 , sin 2 26> 23 = 1 ). 
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(X disappearance x appearance 




10" 1 10° 10 1 10 2 10 3 10 4 10 5 10" 1 10° 10 1 10 2 10 3 10 4 10 5 

UE in eV 2 L/£in eV 2 

Figure 4: The P MM (fig. Qa) and P^ T (fig. ^b) that give the best SK fits. Continuous blu line: standard — > v T fit. 
Dotted red line: — > ^kk fit with intermediate £ = 1/2. Dashed green line: — > i/ r , j/kk /z£ Zarge £. 



of disappearance experiments. For example the on-going K2K experiment [17] should observe only 65% -r 
85% of the events with respect to the no-oscillation case, while the larger range 30% -j- 90% is allowed by a 
standard fit of SK data. Hence K2K could put the constraint R < 0.01 mm on the scenario under study. 

3. A related characteristic feature of the same survival probability is a precocious disappearance of (and 
appearance of v T if V t3 ^ 0) at relatively low LfE v . The effect is not big since otherwise it would have made 
difficult the same SK fit or it could have been in conflict with the null result of CHORUS and NOMAD JTs). 
At small L/E v , an analytic approximation for P a p is 




Note the linear dependence on L, rather than the quadratic one characteristic of the standard oscillation 
formula. The effects related to (|l^) arc not without interest already for the downward going in SK, in 
view of the parameters shown in fig. |^. 

4. A partial suppression of the ^-events might occur in an appearance experiment trying to measure P^ T at 
large L/E u relative to the expectation in the case of the standard — > v r interpretation of the data. 

5. Earth matter effects suppress the oscillations of a long-baseline beam of with large energy E v > lOGeV, 
although with a slower E v dependence than in the usual case. On the contrary earth matter effects enhance 
the oscillations of a beam. 

All these points deserve further quantitative investigation. 



8 Constraints from nucleosynthesis 

So far we have compared the phenomenology of neutrinos from extra dimensions with oscillation experiments. In 
this section we study the constraint arising from standard big-bang nucleosynthesis. As well known, the danger 
for nucleosynthesis is that too many KK neutrino modes are produced before the time of nucleosynthesis. Two 
different production mechanisms have to be studied: (i) by incoherent scattering; (ii) by coherent oscillations. 

In short, the incoherent production is not dangerous provided a temperature T* is assumed, at which the 
abundance of relic KK neutrinos is negligibly small. T* is in the GeV range for 5 — 2 and can easily be much 
higher as 5 increases. This agrees qualitatively with several statements in the literature j|, 

The case with coherent oscillations is more delicate since the evolution of an infinite number of different neutrino 
mass cigenstates is not easy to follow, especially in the case of a sizable £ parameter. 
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A simple way to get an idea of what happens is the following. Since the oscillation frequencies arc large relative 
to the collision or the expansion rates at any temperature T close to nucleosynthesis, one can estimate if any one 
of the mass eigenstates may reach equilibrium by comparing to the expansion rate H{T) the effective interaction 
rate of the n-th state T n (T) w U§ n T u (T). Here T V (T) is the typical interaction rate of a standard neutrino with 
the medium and Uq u is given in (jioj) in terms of the matter potential p. If one neglects any matter asymmetry in 
the primordial medium, either original or generated by the neutrino evolution equations themselves []19| 



p = -2c(E v RG F ) 2 T 4 /a e 



(16) 



where a om and Gf are the fine structure and Fermi constants, E v is the neutrino energy, to be suitably averaged, 
and c is a numerical coefficient close to 0.2 for p and r neutrinos. The p parameter has a negative fixed sign both 
for neutrinos and antineutrinos. As discussed in section |[ a negative p implies a MSW suppression of T n (T). Fig. [s| 
shows the plot of T n (T)/H(T) versus T for £ = 3 for a large number of mass eigenstates different from the special 
one, which coincides with the interacting neutrino at large temperature. The rise at small T is a manifestation 
of the T 3 -behavior of T v /H , whereas the fall at large T is a consequence of the MSW suppression of the effective 
mixing factor [/g n . None of the states is close to equilibrium at any temperature. The opposite would have been 
a clear problem for nucleosynthesis since the typical temperatures are close to the MeV. It is harder to say if the 
situation described in fig. |B| is compatible with standard nucleosynthesis. Even though the energy density stored 
in any state is small because none gets even close to equilibrium, their total energy density might be significant. 

To answer this question requires solving the coupled sys- 
tem of infinite equations for the neutrino density matrix p. 
Here we only check that standard big-bang nucleosynthesis 
is not grossly inconsistent with the picture described above. 
In the T-dependent W ma ttor (T) eigenstatc basis, neutrino 
propagation averages to zero the off-diagonal elements of p 
much faster than the typical interaction or expansion rates 
at T ~ MeV. Setting to zero the off-diagonal elements 
Pnm, ti ^ m, and assuming that neutrinos are in kinetic 
equilibrium, the evolution equations for the neutrino densi- 
ties N n = p nn of all individual neutrino mass eigenstates, 
normalized to the equilibrium density, are 
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0.02 
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dlnT 



r,(r) 

H(T) 



U* n (T)(l-N n ) 



(17) 








We have explicitly included annihilation processes and in- 
dicated with • • ■ additional terms due elastic scatterings. 
These terms have comparable rates and redistribute the to- 
tal neutrino number N in its various components, but do 
not directly affect the total neutrino number N = J^n^n- 

Summing all the equations we find the evolution equation for the total neutrino density 
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Figure 5: Values of T n /H for £ 
2.5 l(T 3 eV for n= 1, . . . 100. 



2 and 1/R 



dN 
d\nT 



H(T) 



(18) 



where N v = J2 n Uo n N n is the number density of the interacting neutrino in the flavour basis. If neutrinos start 
with a normalcy distribution at some temperature T, oscillations rapidly convert it into N n (T) = J7q„(T) so that 



n 4 



P vv . Using this value of N v (T) and the approximation P vv 



p, or large T, we estimate the right-hand side of (|18[) as 



II 4 



dN 



r„(T) < 2 



dlnT H(T) y^T) 



1 — 7r£ / y/—p valid at large 



(19) 



As for the individual rates T n /H, the right-hand side of ( |T9| ) is small at T< 1 MeV and progressively increases. 
However, unlike T„/H . it does not get to a maximum but rather flattens out to a horizontal asymptote at a level 
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Af (£ 2 / R) / (10~ 2 eV). In the present case where v is a combination of and v T , the numerical coefficient J\f is of 
order unity and has a large uncertainty. Therefore we cannot conclude that values of £ 2 /R as 10~ 2 eV, relevant to 
the atmospheric neutrino anomaly, give neutrino oscillations incompatible with standard big-bang nucleosynthesis. 
A detailed and non trivial study could give significant constraints. An unambiguous problem, if present, could not 
be avoided as usual by invoking a large lepton asymmetry which, in the case under study, can only suppress v or 
v oscillations, but not both. 



9 Constraints from supernova 1987a 

In this section we show that the SN1987a observation is likely to put a very severe constraint on i?, not compatible, 
if taken at face value, with the phenomenology of neutrino oscillations as discussed in the previous sections. This 
is because a large number of KK neutrinos would carry away from the supernova (SN) too much energy in invisible 
channels, thus shortening the observed neutrino burst in an unacceptable way. More precisely, we demand that the 
energy loss rate in invisible channels, W- mv , is less than about 10 19 erg/gr • sec for typical average conditions of a 
SN core with a density p as 3 • 10 14 gr/cm 3 and a temperature T as 30MeV J2(J. As in the case of nucleosynthesis, 
one has to consider the production of KK neutrinos both by incoherent scatterings and by coherent oscillations. 

The rate of incoherent production of a single KK state is estimated as rj nc as [rn/E) 2 T iyi where to is the largest of 
the masses rm introduced in section 2, E ss 100 MeV is a typical neutrino energy and Y v as Gpn^E 2 is the collision 
rate of a standard neutrino in the supernova core, in terms of the nucleon density rijv- Accounting for the number 
of KK states, approximately RE, the corresponding energy loss rate can be estimated as Wj nv ~ (m/E) 2 ER ■ W v , 
where W v as 10 27 erg/gr • sec is the energy loss rate produced by T v if the standard neutrinos were not trapped. 
Requiring Wi nv /W u < 10~~ 8 gives (10 -3 eV • R)(m 2 /10~ 3 eV 2 ) < 1 which is not incompatible with the values of m 
and R considered in the previous sections. Note, however, that had we considered 5 > 1 extra dimensions, all with 
comparable radii, the increase in the number of KK states of energy less than E, about (RE) S , would have lead to 
a stronger constraint on R 

R < 0.2 mm • lO" 11 ^ 1 /*) ( ^ , 

VlO- 3 eV 

Let us now consider the effect of coherent oscillations of the active neutrino state associated with the mass to into 
its KK tower. Depending on the sign of /i in the supernova core, MSW resonant oscillations will take place either 
in the neutrino or in the antineutrino channel. Since the neutrino effective potential V in the core is of order 

10 eV, the KK states interested in the resonant transitions have a mass m n ~ n/R as {VE) 1 / 2 as 10 4 ^ 5 eV, with a 
number of states interested in the resonant transition of order n ros as R ■ 10 A ^ 5 eV. The effect of these oscillations 
can be reliably estimated for small £ = mR. This is so because in such a case, due to the small vacuum mixing 
angle between the n-th KK state and the standard neutrino, 9 n as £/n, the width of the n-th resonance is smaller 
than the separation between two contiguous resonances. As a consequence, the survival probability for a standard 
neutrino (or antineutrino) produced in the core can be computed as the product of the survival probabilities in 
each resonance, P vv as Y\ n P n @- Furthermore, for each P n one can take an approximate formula describing the 
transition with initial and final densities respectively much larger and much smaller than the resonance density, 
P n as e"^™/ 2 , where 

4£ 2 V 
7,1 ~ R?E dV/dr rcs 

is the adiabaticity parameter at the n-th resonance crossing. Approximating V/(dV/dr) with the radius of the 
core, r coro as 10km, j n is of order m 2 /10~ 3 eV 2 and is approximately n-independent. Hence the individual P n are 
expected to deviate sensibly from unity, which is a fortiori the case for the overall probability P„„ due to the large 
number of resonances crossed. This is untenable no matter what the fate will be of the KK states produced in the 
oscillations, either escaping or trapped in the core, since in both cases the luminosity in standard neutrinos will 
be drastically reduced. We do not even see how such a conclusion could be avoided by having 1/R approaching 
to, or £ > 1, even though a precise calculation in this case is not easy to make due to the overlap of the different 
resonances. 

We therefore believe that this sets an upper limit on R close to l/10 4 ~ 5 eV. Above this value no resonance is 
crossed anymore and the mixing angle of any KK state with the standard neutrinos becomes negligibly small. This 
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bound, together with (|2Cj), would seem to indicate that independently from the number of dimensions, the largest 
of the radii cannot be bigger than about 1A. 



10 Conclusions 

In the discussion about possible patterns of physics beyond the SM, the existence is being considered of some large 
extra dimensions where SM singlets could propagate. The possibility of lowering the scale where gravity becomes 
strong is a striking consequence of this hypothesis. Hence, the hierarchy problem is generally quoted as its main 
motivation. On the other hand, unlike the case, e.g., of supersymmetric unification, we know of no observation, 
direct or indirect, that supports this picture. The problem of finding possible signals of large extra dimensions is, 
therefore, particularly acute. 

In this paper we have shown that neutrinos offer a concrete possibility of testing the existence of large extra 
dimensions, alternative to those considered in high energy collisions or in gravity measurements. This is so if there 
is a large extra-dimension that extends above 0.1/im, where right-handed neutrinos propagate. We think that 
neutrino oscillation experiments can fully explore this range of distances. Constraints arise in principle on this 
phenomenology both from standard big-bang nucleosynthesis and from SN1987a. As we have discussed, particularly 
constraining appears to be the consideration of the neutrino luminosity from SN1987a. If taken at face value, a 
maximum radius of any extra dimension is implied by these considerations close to lA, making impossible any 
observation in neutrino oscillation phenomena. This bound is likely to be valid, even though we cannot perform 
a clear calculation in all the relevant range of the parameters. However ad hoc ways to evade it without affecting 
the signals discussed in sections 6, 7 may exist, like suitable decays of the heavier KK neutrinos in the 100 KeV 
range, which undergo the resonant transitions in the supernova core. In our view, it is anyhow worth pursuing a 
direct comparison with neutrino oscillation data. 
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A Explicit diagonalization of the Hamiltonian 

From (||) we arrive to the Dirac mass matrix (Q) with the following definitions: 



for n > 1 



and 



(21) 



In this way ipf and the linear combinations orthogonal to i/^ and i/j? decouple and can be forgotten. The resulting 
effective Hamiltonian (^) for the first fc max KK states is H c g = X/2E V R 2 where the dimensionless matrix X has 
an 'almost diagonal' form: keeping only the first fc max KK excitations its only non-vanishing off-diagonal elements 
are X k0 = X ok = k\/2^ (k > 1). Its diagonal elements are Xkk = fc 2 for fc > 1 and Aoo = (2fc max + 1)£ 2 + fj,, where 
the term parameterizes matter effects. Calling A 2 the eigenvalues of X, the eigenvalue equation can be written 
in the closed form 



detpf - A 2 1] = 



(A 00 -A 2 )-£ 



XkoXok 



fc=i 



X kk - A 2 



IK* 

k=l 



In the limit fc max -^oowe obtain the eigenvalue equation ( 
matrix may have negative eigenvalues: for [i < £ 2 , in fact, 



A 



A 2 



A 2 



fe=i 



A 2 



n ( fc2 - a2 ) 



fc=i 



). As already remarked, when fj, < the squared mass 
t) admits one A 2 < with A purely imaginary. 



For any eigenvectors e^(A) the secular equation gives a relation between the various elements and the first one: 



e k/ e o = — fc\/2£/(fc 2 — A 2 ). Using the normalization condition 



1 wc arrive at 



T7T- = = 1 + 2C 



U 2 



4(K) 



2 \ - ft, 

^(fc 2 -A 2 ) 2 



2 + 7T 2 £ 2 (l+COt 2 7rA„) - 



n 2 e 



cot 7rA„ 



(22) 



which, using (p|), gives (|T(]). It is not difficult to obtain simple expressions that extend the large-£ oscillation 
amplitudes (p"lj) and ( ]l2| ) to include matter effects. 
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If p < 0, one has to separately include the contribution of the special cigenstate with A 2 < 0. When p <C — (7r£ 2 ) 2 
this state has a small KK component and P VU {L — ► oo) « C/q w 1 — ?r£ /y/— \i so that matter effects suppress 
oscillations. On the contrary when p > cq. ( |l2] ) remains valid. Moreover when the matter potential is very large, 
p = /i/(7r£ 2 ) 2 3> 1, the effective L/E v above which neutrino oscillations are significant is reduced by a factor p 1 / 2 
with respect to the p = value given by eq. (O). 
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